#!/usr/bin/env python
# -*- coding: utf-8 -*-
"""
省份地理边界数据提取主脚本

此脚本是对gis模块功能的简单封装，保持与原始脚本相同的功能，
但使用了模块化的设计，便于维护和扩展。

使用方法与原始脚本相同，但内部实现更加结构化。
"""

from scripts.gis.province_utils import process_province_data


def main():
    # 指定要提取的省份列表
    target_provinces = [
        '黑龙江', '新疆', '山西', '宁夏', '山东', '河南', '吉林', '辽宁', 
        '天津', '青海', '陕西', '内蒙古', '河北', '北京', '甘肃'
    ]
    
    # 读取省份shapefile文件
    shp_path = r'd:\Workspace\Meicao\wangrui\dry-wet-abrupt\GS(2023)2767审图号\省面.shp'
    
    # 设置输出目录
    output_dir = r'/'
    
    # 处理省份数据
    result = process_province_data(
        shp_path=shp_path,
        target_provinces=target_provinces,
        output_dir=output_dir,
        csv_filename='province_bounds_detailed.csv',
        map_filename='province_map.png',
        html_filename='province_report.html'
    )
    
    print("\n输出文件:")
    for key, path in result.items():
        print(f"- {key}: {path}")


if __name__ == '__main__':
    main()